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I. INTRODUCTION 

The increasing experimental capabilities in the fabri- 
cation of semiconductor nano-crystallites (quantum dots, 
QD) for a wide spectra of sizes drives the experimental 
and theoretical interest for the overall comprehension of 
these systems since, due to the quantum confinement, the 
electronic and optical properties are strongly dot size- 
dependent. El In particular, it is well established that the 
optical band-gap energy is augmented (as compared with 
bulk values) towards higher energies when the dot size 
decreases. This behavior is the result of two opposite ef- 
fects: while the single-particle band gap energy is shifted 
to higher energies, the binding energy of the electron- 
hole pair (exciton) created by photoexcitation adds an 
attractive correction which gives rise to a net blue shift 
but with a weaker dot size dependence as compared with 
the corresponding to the single-particle band gap energy. 

The first and simplest theoretical approach used in the 
analysis of electronic properties of nanaatructures was the 
effective mass approximation (EMA),0i3 which in princi- 
ple is capable of covering the full range from large to small 
quantum dots, with the length scale given roughly by a 
few lattice constants. It is then clear that there should 
exist a critical dot size for the validity of the EMA, below 
which the approximation (essentially based in the semi- 
conductor bulk properties) must fail. Fortunately, in this 
limit of very small crystallites (up to a maximum diam- 
eter of about 30 A) other more microscopic n a [ Dproaches 
such as the empirical tight-binding (F^B^,mD the em- 
pirical pseudopoteptial method (EPM)j3~1i3 and ab initio 
pseudo-potentialEU calculations are possible and improve- 
ments in quantitative agreement with experimental re- 
sults has been reached in this small size regime. Based in 
a comparison between their results for the single-particle 
band gap energy, the Coulomb and exchange excitonic 
corrections and the corresponding EMA results, most of 
the authors of this group conclude in disregard the EMA 



as a quantitative, and even qualitative method for the 
study of nanocrystallites. However, we must put atten- 
tion to the fact that in most of the cases the EMA re- 
sults were obtained under the assumption of perfect con- 
finement (infinite barriers, IEMA) for both electron and 
hole. Several authors pointed out that this extreme hard- 
wall constraint (which is not inherent to this approach) 
is the main reason for the discrepancy with more micro- 
scopic approaches and carried out calculations assuming 
an incomplete confine men t (finite barriers, FEMA) for 
the electron and hole.E3Ej The output of such calcula- 
tions is that while the single-particle band gap energy 
is not well described by FEMA, excellent quantitative 
agreement was obtained for the excitonic properties with 
the more sophisticated calculations available at present . 
This can be understood from the fact that while single- 
particle properties depend directly on the fine structure 
of the electron and hole envelope functions, exciton prop- 
erties are given by averaged (integrated) expectation val- 
ues with these envelope functions (see below). As a con- 
sequence, excitonic related properties are much less sensi- 
tive to fine details of the electron and hole wave functions 
than single particle properties. The encouraging results 
obtained in this way allow us to keep FEMA as a quali- 
tative, flexible and versatile theoretical tool for the study 
of the excitonic physics of quantum dots. 

Another question to analyze in detail in order to be 
able to make a proper comparison between theory and 
experiment is the change in the excitonic energy due to 
the dielectric mismatch between the quantum dot and the 
surrounding matrix. This brings immediately the issue 
that electrostatic screening in semiconductor quantum 
dots is expected to be different from its bulk counter- 
part. In particular, reasonable doubts arises about if a 
macroscopic treatment based on the use of an effective 
dot dielectric constant (e) is still meaningful. This topic 
was recently addressed-in three different calculations, a 
phenomcnolpgical onejlj and two full quantum mechan- 
ical onesffij All of them conclude that the confinement 
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should reduce e, although differed with respect to the 
magnitude of this reduction. Besides, and this is the key 
point, in the last two microscopic calculations, through 
a careful comparison between microscopic and macro- 
scopic results for the electrostatic screening, the validity 
of macroscopic type treatments was established, although 
with suitable generalizations. This size-dependent effect 
on the dot dielectric response can be easily incorporated 
to our approach, but we shown at the end of Section III 
that this correction is quantitatively quite small as com- 
pared with the size-dependence effects on the Coulomb 
and self-energies. Consequently, in the rest of the paper 
we ascribe to the dot a bulk dielectric constant, which 
is an excellent approximation for not too small quantum 
dots (dot diameter, d > 20 A). For simplicity, most of the 
times the effect of the dielectric mismatch is incorporated 
assuming that the dielectric constant changes discontin- 
uously from its value in the quantum dot to the value 
corresponding to the matrix (step-like model, see Fig. 
la). In this case, and according to the classical image- 
charge concepts of basic electrostatics,E£l the presence of a 
charged particle inside the quantum dot induces charges 
at the boundaries which have the same sign as the source 
charge if the dielectric constant inside the dot is higher 
than outside (the general situation). Within IEMA, this 
gives rise to a repulsive (positive) potential energy corre- 
sponding to the particle interacting with its own induced 
charge; this is usually denoted as a self-polarization en- 
ergy correction. Besides, when we study excitons, the 
polarization effects also modify the direct Coulomb in- 
teraction since the induced charge by the electron (hole) 
can attractively interacts with the hole (electron) for the 
usual case. Consequently, there are two opposite contri- 
butions to the excitonic energy leaded by polarization 
effects. For spherical quantum dots, and assuming i) 
infinite confinement for the particles (IEMA) and, ii) a 
step-like model for the dielectric profile, it is possible to 
demonstrate' that both effects are almost exactly com- 
pensated and then, the exitonic energy is reduced to the 
bare electron-hole attraction. In a previous work, we 
have verified tiiat the same cancelation occurs for cubic 
quantum dots.cj However, as soon as we depart from this 
two strong assumptions, the compensation between the 
polarization effects will not verify anymore. Therefore, 
the two polarization contributions to the energy must be 
taked into account in more realistic model calculations. 

A direct consequence of the step-like dielectric profile 
is that all the induced charge is perfectly localized at the 
surface boundary, leading to a divergent self-polarization 
potential at the dot surface. This feature is not compat- 
ible with the electron and hole incomplete confinement 
associated with finite potential barriers, since in this case 
both the particle and its induced charge could coincide 
at the dot bojAtidary and then, the self-polarization en- 
ergy diverges.tll This unphysical result clearly show that 
the standard electrostatic description of the interface be- 
tween dielectric media fails at distances comparable to 
the interatomic distance. The difficulty is automatically 



solved in the infinite barrier case: as a consequence of 
the hard-wall boundary condition the electron and hole 
envelope wave functions are zero precisely where the self- 
induced potential diverges (the dot boundary), and then 
the self-polarization energy remains finite. For the fi- 
nite barrier case, the anomalous behavior can be "regu- 
larized" by introducing a phenomenological cut-off dis- 
tance of ±he-&ame order of magnitude as the interatomic 
distance£a~E3 However, this proposal of modifying the 
functional form of the potential does not follow from 
any calculation; it is only a mathematical trade-off de- 
signed to avoid the divergence. .-A more rigorous ap- 
proach was carried out by SternEil who, in the study 
of a semiconductor-insulator Si/Si02 planar interface, 
solves the problem proposing a model dielectric inter- 
face in which the dielectric constant changes continuously 
between its extreme values (see Fig. lb for our three- 
dimensional analogue). In this way the induced charge 
is spread along the dielectric interface region and the di- 
vergence in the self-polarization energy disappears, as it 
physically should be. The consideration of an interface 
with finite thickness is in parallel with the situation in 
real samples because the experimental growth procedures 
do not guarantee a sharp profile without interdiffusion 
between dot and matrix materials. 



Figure 1 P.G. Bolcatto and C.R. Proetto 





(a) 






R r 

(b) 

'\ 

NS. 

— ^ 


R 


-5 R + S r 
(c) 







R-5=R R,_, R, R N _, = R + S ' 



FIG. 1. Spatial profiles for the dielectric interface, (a) 
Step-like model; (b) linear (dot-dashed line) and cosine-like 
(full line) models; (c) schematic discretization by N— steps 
of an arbitrary model of the dielectric interface. The step-like 
model corresponds to the choice 8 = 0, N = 1, while in the 
limit N — > oo the profile of Fig. lc can be made arbitrarely 
close to any of the continous profiles of Fig. lb. 
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As is discussed above, the evaluation of the optical 
band gap requires two independent calculations: the 
single-particle band gap energy (blue shift) and the exci- 
tonic corrections (red shift). By using EMA, the accurate 
calculation of the former contribution requires considera- 
tion of other ingredients as a multiband-k»p finite barrier 
calculation of the dot band structure,E2l which is beyond 
the scope of the present work. Therefore, the aim of the 
present paper is to give insights about the excitonic cor- 
rection in quantum dots when the combined effects of 
partial confinement by finite barriers and dielectric mis- 
match are considered. 

The rest of the paper is organized as follows: in Sec- 
tion II we give the basic formulae used for the calcula- 
tion of the excitonic energies, in Section III we discuss 
the results, and finally in Section IV we summarize our 
conclusions. 



II. THEORY 

In a previous work@ we have demonstrated that the 
excitonic properties (even in the presence of large di- 
electric mismatch) are essentially independent of the dot 
shape, so we will concentrate our analysis in the simplest 
geometry: spherical quantum dots with radius R. We at- 
tach to the dot-acting semiconductor a dielectric constant 
Ein, while the outside media where the dot is embedded 
has a dielectric constant e ou t , and we will study the effect 
of different models for the dielectric interface connecting 
both limits (see Fig. 1). Within the envelope wavaiimc- 
tion approach to the effective mass approximation,!^!] the 
relevant Hamiltonian for an electron-hole system is given 
by 

H (r e ,r h ) = H e (r e ) + H h (r h ) +V c (r e ,r h ) , (1) 

where the first and second terms on the right-hand side 
are the single-particle contributions for the electron and 
hole respectively, and the last one is the generalized 
electron-hole Coulomb interaction (to be discussed be- 
low). The single-particle Hamiltonian Hi (r^) (i = e,h) 
is defined as 

^( r O = -^ArrV?+K(r i ) + K)i(r i ), (2) 

where the first term is the kinetic energy for a parti- 
cle with effective mass m., (r,), V s (ri) corresponds to 
the above discussed self-polarization energy that can be 
obtained as a particular case of y c (r e ,r/j) (see below), 
and the last one is the confining potential which we as- 
sume as a spherical quantum well-like potential defined 
by Vbi (rj) = 9 (r-j — R) Voi, being the step function 
and Voi the barrier height. The dependence of mi on n 
arises from the fact that the particles have different effec- 
tive masses depending on their location, inside or outside 
the dot. 



A. Coulomb and self-polarization potentials 

The generalized Coulomb potential between two par- 
ticles with coordinates r and r', and charges e and e' 
respectively, can be calculated through V c (r,r') = e 
3>(r, r'), where $(r,r') is the electrostatic potential 
which verifies the Poisson equation 

V r • e (r) V r $ (r, r') = -Ane'S (r - r') , (3) 

where e(r) is the (local) spatially dependent dielectric 
function. This equation is quite simple to solve analyti- 
cally if we assume that e(r) depends only on r, and also 
that the spatial dependence for e(r) corresponds to the 
step-like model, as shown in Fig. la (for instance, see 
Refs. lb and 3). The important points in the derivation 
of such solution are: first, under the assumption of spher- 
ical symmetry Eq.(3) is separable in angular and radial 
parts, being only the resulting radial differential equation 
the non-trivial part to solve; second, the solution of this 
radial equation is quite easy to find if e(r) is constant. 
For this particular case, the radial solution of Eq. (0) 
consists of product of powers of the coordinates r and 
r'; the step-like model of Fig. la allows then to propose 
for each region at both sides of the interface two differ- 
ent solutions of this type, and finally obtain the solution 
for the whole space through the standard electrostatic 
boundary conditions which $ should satisfy at the inter- 
face. It is also possible to obtain an analytical solution of 
Eq. (|) in the case where e (r) depends linearly on only 
one coordinate, we say z, and is independent of x and 
y, as for example in_thc case of a simple semiconductor- 
insulator junction.EU The disappointing feature of this 
exact solution is that the electrostatic potential exhibits 
a logarithmic divergence at the edges of the transition 
layer, which can be traced back to the discontinuous first 
derivative at these points of the linear model for the di- 
electric interface. The three-dimensional solution for the 
linear model of the dielectric interface, even assuming 
spherical symmetry, is quite complicated to find since 
the Eq. (||) becomes a second-order differential equation 
with spatially dependent coefficients. We will not pursue 
the analytical resolution of this problem, which even if 
feasible, is restricted to just a particular model for the 
dielectric interface. 

Instead, by exploiting the exact solution of the single- 
step case, and following the spirit of Ref. [ EHl, we will 
consider an arbitrary smooth radial profile for e (r) and 
proceed as follows: first we define a finite interface of 
width 25 centered at R. For r < R — S (well inside the 
dot) the dielectric constant takes the value correspond- 
ing to the bulk dot-acting semiconductor, e (r) = £j„. For 
r > R + 5 (well outside the dot) e (r) is fixed at the value 
corresponding to the surrounding matrix, s (r) = £ ut- 
Between them, for R — S < r < R + 5, we can choose 
any analytical and physically plausible function that con- 
nects these extreme values as is pictured in Fig. lb. Af- 
terwards we subdivide the dielectric interface in N — 1 
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regions and in each one of them we approximate the di- 
electric function by a constant value (typically the mean 
value of the selected function in this zone). Thus, the 
original continuous profile is approximated by an A-step 
one in which each single-step problem has a known so- 
lution. This approach allows us to connect £,„ and e ou t 
with any arbitrary function being the number of steps A 
the only parameter to adjust, the condition being that 
it should provide a good approximation to the true con- 
necting function. Strictly, our approach is such that when 
A — > oo we recover the exact solution of Eq. (3). As a re- 
sult of this procedure, the expression for the generalized 
Coulomb potential is (see Fig. lc to clarify the notation) 



V*> m (r,r') = — 5>(cos 7 ) 



F im ,i (r, r') 
(1 - Pm,lQm,l)' 



(4) 



where the superindex i ( m ) indicates that the parti- 
cle with coordinate r ( r' ) is placed at the i— region ( 
to— region); Pi (cos 7) are the Legendre polynomials of or- 
der I and 7 the angle between both particles (measured 
from the origin at the dot center). With i and m taking 
values between and A, the Coulomb potential V c (r, r') 
is defined in the whole space. The functions Fi m j (r, r') 
are given by 



Fimj (r, r') = [p n 



(i+i) 



m-l 1 D-(2i+X) 



1 + qj,lR. 



-(21+1) 



(5) 



if % < m, 



F imj i (r, r') = p m> 



-(l+l) 



r< + q m ,ir. 



(6) 



for i = to and 
Fi mt i (r, r') = ^pijr 



-('+!) 



"< + q m ,ir < 



» 1 -L n r?2i+l 

J=m +i 1 +Pj,i R j-i 



(7) 



when i > to. In the previous equations r < (r>) is the 
smallest (greatest) distance between r and r'. The gy 
coefficients are defined so that qo i — (VZ) and for 
< i < to are expressed through the following recur- 
sive expressions 



qu 



J 21 ' J 22 
rti,l 1 rii.l ' 

*11 + 3i-l,i<-'l2 



(8) 



where the S— matrix is 
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l(e 



(/ + 1) + £i-lZ 

\ p2i+l 
•i-lj 



(/ + 1) (£< - ft 
£*Z + £i_! (Z + 1) 



-(2i+l) 



(9) 



Similarly for the py coefficients, pjy,i 
to < i < A 
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Pi+IA&U + fog 
Pi+U^21 + "^22 



(VZ) and for 
(10) 



where the E— matrix is 



£l (/ + 1) + £ i+1 Z (/ + 1) (e, - £ i+1 ) i?-( 2i+1 > 
LZ(e,-e l+ i)i?f +1 £ i+ i (Z + 1) + Ei l 
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Some interesting features of Eqs. (4)-(ll) are: i) the 
coefficients gy and p m ,l are all identically zero in the 
absence of dielectric mismatch at the dot boundary 
(Ei n = £ ou t)\ h) the function F im j(r,r') defined in Eqs. 
(5)-(7) has four different contributions, corresponding 
to particle-particle interactions (the term rL / r 1 ^ 1 ), 
particle-induced charge interactions (the two terms ei- 
ther proportional to or p m ,i), and the induced charge 
- induced charge interaction (the term proportional to 
the product <jy p m ,l)- Only the first one, which rep- 
resents the bare Coulomb interaction, survives in the 
homogeneous limit £j n = e out ] iii) in the calculation 
of V^ m (r,r'), only the qij's (p m /s) with < i < to 
(to < i < A) are necessary. For the exciton problem, 
e' = - e. 

From Eqs. ^ and (||), the self-polarization potential 
can be calculated by taking r = r', e = e', eliminating 
from (||) the direct contribution r l K / r 1 ^ 1 , and dividing 
by 2 as corresponds to a self-energy. The result is 



V s m (r) 



\ " 



2e m j^t (1 - p m ,iq m ,i) 

p m ,ir 21 + Pm,l<lm,ir~ l + q m jr~ 2{l+1) ] . (12) 



With to running from to AT, the self-polarization po- 
tential 14 (r) is defined in the whole space. 

It is straightforward to check that in the limit of S = 
(A = 1) Eqs. (g) and @ reduced to the usual expres- 
sions for and V s used in the studKXif spherical QD's 
with boundary dielectric mismatch.o&Ej 

The choice of different profiles for the dielectric inter- 
face modify both the strength and the functional form of 
the potentials. In Fig. 2 we show the self-polarization 



potential (Fig. 2a) and the 



fcr 



component of the 



generalized Coulomb potential (Fig. 2b) for a small 
quantum dot where R = 10 A, £j n = 12.6 (GaAs), 
£out — 1 (vacuum). The dot-dashed curves correspond 
to the case of the step-like model for e (r) (that is S = 
and N = 1). The self-polarization potential for this case 
exhibits the already discussed divergent behavior at the 
interface, which leads to a divergent self-energy for the 
finite barrier situation. The dotted curve is the result 
of choosing a linear model for the dielectric interface be- 
tween ei n and e ou t, and the solid line corresponds to a 
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profile proportional to a cosine function; in both cases, 
8 = 5 A and N = 500. 

Figure 2 P.G. Bolcatto and C.R. Proetto 
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FIG. 2. Self-polarization potential (upper pannel), and 
1 = component of the generalized Coulomb potential (lower 
pannel) as a function of r, for a quantum dot with R = 10 
A, and several models of the dielectric interface. Dot-dashed 
line, 5 = 0; dotted line, S = 5 A (linear profile); solid line, 
8 = 5 A (cosine- like profile). For the generalized Coulomb 
potential, the source charge is located at r' = 2.5 A. 

We can observe in Fig. 2a that the consideration of a 
finite interface [independently of the form of e(r)] elim- 
inates the divergence present when 8 = 0. It is worth 
to emphasize here the important difference between the 
"regularized" self-polarization potential and the solutions 
which arise from the finite size dielectric interface. Faced 
with the dilemma of a divergent self-polarization poten- 
tial and accordingly a divergent self-energy, no physical 
guidance (except that the final result should be finite) 
helps one in choosing a regularization procedure for the 
step-like model of the dielectric interface. In our ap- 
proach, we do have such a physical guidance: the di- 
electric function must change smoothly in moving from 
the dot towards the matrix (presumably on the scale of 
a few A), and from this physically plausible assumption 
the non-divergent self-polarization potential emerge nat- 



urally from calculation. The lack of a physical guidance 
for the case of a step-like model of the dielectric interface 
is clearly reflected in the fact that the "regularized" self- 
polarization potentials of previous worksEj~tj are quite 
different from our physically correct self-polarization po- 
tentials shown in Fig. 2. For instance, a regularized 
potential following the criteria of Ref. [ should be 
made by replacing the dot-dashed curve for 5 A< R < 
15 Aby a straight line connecting the black dots. 

Similarly to the planar semiconductor-insulator inter- 
face, the linear model for the dielectric interface leads to 
a potential which exhibits singular points at the left and 
right edge of the interface .ED precisely where the deriva- 
tive of s(r) is not continuous. While from our numerical 
approach it is difficult to determine the type of singu- 
larity, it can be integrated easily numerically, so we can 
safely state that the singularity is integrable. However, 
a cosine-like function seems to be a more physical choice 
because it gives rise to a well-defined self-polarization 
potential everywhere. An important feature of this po- 
tential for all cases is that it can be either repulsive or 
attractive. When the source charge is located in the re- 
gion with greater dielectric constant, the induced charge 
has the same sign as the source and then the interac- 
tion is repulsive. On the contrary, if the source charge 
is placed in the zone with lower dielectric constant, the 
induced charge has opposite sign and the interaction is 
attractive. Note however the strong asymmetry between 
the repulsive and attractive potential strengths. In next 
section, we will see that this singular behavior has a rel- 
evant role on the calculation of the excitonic energies. 

Regarding the generalized Coulomb potential V c , 
all the curves have the universal expected r _1 decay 
screened by e out outside the dot (r > R + 8) . The results 
for linear- or cosine-like profiles for e (r) are quite similar 
(dotted and solid lines, respectively), being the potential 
strength in both cases smaller as compared with the step- 
like model. This can be explained as follows; from Eq. 
(4), the I = component of the generalized Coulomb po- 
tential for the step-like model for the dielectric interface 
is 



V c °>°(r,r') = 



R 



^out 



(13) 



for r < r' (the region where the potential is constant 
of Fig. 2b). The first term in Eq. (13) corresponds to 
the bare Coulomb potential, the second to the potential 
generated by the charges induced at the dot boundary. 
Note that this term is proportional to the difference of 
dielectric constants at both sides of the interface; conse- 
quently, by giving a finite size to the dielectric interface 
this term can only decrease in magnitude, as the dif- 
ference between two successive dielectric constants ap- 
proach zero for N — * oo. Already from the results pre- 
sented in Fig. 2b a decrease of the Coulomb energy of 
the exciton for the case of a finite size dielectric interface 
can be anticipated, as compared with the 8 = interface. 
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B. Zero-order eigenvalues and eigenfunctions 

As we will consider quantum dots with characteris- 
tic dimensions smaller than the exciton effective Bohr 
radius, we will use the strong confinement approxima- 
tion (SCA) for the calculation of the eigenvalues, wave 
functions, and excitonic energies .1313 In this approach 
the self-polarization and Coulomb interactions (which 
scales as are considered as a perturbation against 

the kinetic energy (which scales as i?~ 2 ). Accord- 
ing to this, we propose a separable wave function for 
the electron- hole system ^ (r e ,Th) — ip e (re) iph ( r h) = 
Rni (r e ) Y lm (6 e , tp B ) Rnn> (r h ) Y Vm , (9 h , tp h ), where 
Rniifi) are the radial component of the wave-function, 
and Yi m {9, ip) are spherical harmonics. As a consequence, 
the resolution of the zero-order Hamiltonian consists of 
two equivalent problems 



H<? ] (r,)^ (r 4 ) 



2m, (n) 

Eitpi (Ti) , 



V 2 + V 0l (r,) 



i>% (Ti) 

(14) 



where i = e,h. 

The complete solutions of Eq. (14) turned out to 
be spherical Bessel functions of real and imaginary 
arguments ;E3 focusing the attention on the ground state 
(n = 1, I — m = 0), the angular contribution of ipi (r<) 
becomes a constant and the normalized solution of Eq. 
(0)is 



i>% (t^ 



A, 



&(R-n) 



sin (k, 



in , 



out„ 



e (n - R) b, 



exp (-k° ut n) 



(15) 



where the eigenvalues are defined through k\ n = 
l E l lh 2 fl 2 and k° ut = [2m? ui (V 0i - /fe2]1 2 



(2m\ 



/K 1 



As the electron (or hole) has different effective mass 
depending on the material, the wave function jnust 
verify the Ben Daniel-Duke boundary conditions,!!!! i.e. 
i/h(Bt) = ipi(R + ) and ^(R^/mf 1 = i>[ (R+) /mf*. 
By applying this special boundary conditions we arrive 
to the following implicit eigenvalue equation for the quan- 
tum dot energies 



kTRcot ( k) n R) = 1 



rn) 



(16) 



Ai and Bi are constants determined by normalization 
requirements and are equal to A; = 



-1/2 



and 



[R/2 - sin (2kj n R) /4fcf + sin 2 (kj n R) /^fcf 
Bi = exp (k° ut R) sin (fc™i?) . An important point in this 
treatment is that we maintain a sharp-profile (step-like 
model) for the electronic interface, choosing its bound- 
ary at the midpoint of the dielectric interface; a more 
rigorous treatment should also give a finite size to the 



electronic interface. Our approximation is justified for 
not too small QD's, as then the electronic interface is 
a small fraction of the total dot size, and accordingly 
zero-order eigenvalues and wave-functions are weakly de- 
pendent on the finite size of the interface (typically a few 
angstroms). 
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FIG. 3. Bound-state eigenvalues (thin lines) of the quan- 
tum dot spherical well, for the step- and cosine-like models 
(thick lines) of the electronic interface. For S — (8 7^ 0) 
the eigenvalues are exact (approximated) , as explained in the 
text. For this particular dot size, only one bound state exist 
in the spherical well. 

In order to give a quantitative support to this approxi- 
mation for the electronic interface, we present in Fig. 3 a 
calculation of the first-order correction to the square- well 
eigenvalues, when the confining potential is modified in 
correspondence with the finite size of the dielectric inter- 
face. In particular, we adopt for the finite-size electronic 
interface the same cosine-like profile as for the finite-size 
dielectric interface. Accordingly, full, dashed-dotted, and 
dotted lines correspond to the cosine-like electronic in- 
terface with 8 = 0, 2, and 5 A. In the last two cases, 
the eigenvalues were obtained treating as a perturbation 
the difference between the cosine- and step- like electronic 
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profiles. In the worse case of 5 = 5 A, the correction 
amounts to about 10 %, while for 5 — 2 A (not shown) 
is much smaller (~ 3 %). Similar differences could be 
expected for the perturbed wave functions. Is interest- 
ing to note that as the perturbed eigenvalue is above 
than the corresponding to the step-like electronic profile, 
one can also expect a greater penetration in the barrier. 
However, the decay of the perturbed wave function starts 
before than the decay of the unperturbed one, and con- 
sequently some compensation exist between both effects. 
Finally, the example displayed in Fig. 3, with 5/R = 0.5 
(for 5 = 5 A) can be considered as the less favorable con- 
figuration for this approximation, as the error decreases 
with 5/R. For instance, for 5/R = 5 A / 20 A = 0.25, the 
difference between perturbed and unperturbed eigenval- 
ues is only about 5 %. 

C. Exciton Coulomb energy 

The exciton Coulomb energy is defined as the correc- 
tion to the single-particle size-dependent band gap en- 
ergy needed to create an electron-hole pair inside the 
quantum dot. Then, this excitonic energy includes the 
generalized Coulomb interaction between the negative 
and positive charges as well as the self-polarization ener- 
gies corresponding to the electron and the hole, 

E ex = E e + £/j — Ecoui- (17) 

The Coulomb energy Ecoui and the self-polarization ones 
E e and E^ are calculated within the SCA by performing 
the average value of V c and V s [Eqs. (^) and (|l^), re- 
spectively] using the zero-order wave functions given by 

S< = J dr t i,* (r,) V, (n) ^ (r*) , (18) 

Ecoui = - J dr e dr h \ i\) e (r e )| 2 V c (r e ,r h ) \ip h (r h )\ 2 . 

(19) 

Defined in this way, Ecoui for the electron-hole pair is 
always positive. For the finite barrier case, Ej can be 
either positive or negative, while for the infinite barrier 
case Ej is always positive. 

TABLE I. Material parameters used in the calculation. 
m/mo e Vb (eV) 

in out in out 
electron 0.067 1 12.6 1 4 
hole 0.13 1 12.6 1 oo 



III. RESULTS 

The results of this section were obtained using the ma- 
terial parameters summarized in Table I, corresponding 
to a GaAs QD immersed in vacuum. We choose this sys- 
tem by several reasons: i) the choice of vacuum as the 
surrounding media simplifies considerably the material 
parameter choice, as then Vbe is given by the electron 
affinity of the semiconductor dot material in bulk, while 
for our simple two-band semiconductor model, Voh — ► oo; 
ii) the high dielectric mismatch between GaAs and vac- 
uum maximize the effect of the induced charge at the dot 
boundary in which we are interested. An important point 
to remark is that all the constants characterizing the ma- 
terials were extracted from the bulk values and conse- 
quently they are not adjustable parameters. Also, all the 
results to be presented below for the self-polarization, 
Coulomb, and excitonic energies were obtained by using 
the cosine-like model for the dielectric interface; quite 
similar results were obtained for the linear model of the 
dielectric interface. 
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FIG. 4. Exciton self-polarization energies (left pannel) and 
Coulomb energies (right pannel) as a function of dot size, for 
5 = 0.05, 0.5, 2, and 5 A. The inset of Fig. 4b corresponds to 
Ecoui{d) as a function of d" 1 . The full thick line corresponds 
to the infinite barrier limit for electrons (and 5 — 0), and the 
remaining curves are for Vbe = 4 eV. Full circles and open 
squares correspond to the ground-state binding energy of an 
on-center donor impurity. 

The lower panel of Fig. 4 shows the results for the 
Coulomb energy [Eq. (|l~9|)], for several sizes of the 
dielectric interface, as a function of the dot diameter 
d = 2R. As a reference we include a_thick curve corre- 
sponding to the usually quoted resulifl Ecoui{d) — 3.786 
e 2 /dEi n which results from the combination of infinite 
barriers (Vbe ~~ > oo) with step-like model for the dielec- 
tric interface. We can observe that, as was expected 
from the behavior of I — component of V c shown in 
Fig. 2b, this contribution to the exciton energy decreases 
monotonously when the thickness of the dielectric inter- 
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face increases. For a better visualization of the scaling 
with dot size we plot in the inset the variation of Ecoui 
with the reciprocal diameter 1/d. The thin line corre- 
sponds to 5 = 0.05 A which can be taken as a very good 
approximation to the 6 = case. As it was established in 
Ref . [ , as soon as we go apart from the infinite barrier 
approximation, the widely criticized d _1 scaling of IEMA 
is lost. Here, we have demonstrated that it is possible to 
modify even more the size dependence if we consider a 
finite size dielectric interface connecting the QD mate- 
rial and the surroundings. For instance, for S = 5 A the 
scaling of Ecoui with dot size is proportional to d~ Q S6 . 

The exciton sclf-polarization energy given by E e + E/j 
plotted on the upper panel of Fig. 4 shows an interest- 
ing behavior with a non-monotonic dependence on the 
dielectric interface size. As in the lower panel, the full 
thick line corresponds to the infinite barrier limit com- 
bined with the step-like model of the dielectric inter- 
face; accordingly, the difference between the full thick 
line (6 = 0, Vb e — > oo) and the full thin line (S = 0.05 
A, Voe = 4 eV) is due to the effect of the finite barrier 
size for the latter. In particular, polarization effects are 
strongly suppressed for near zero-width interfaces, after- 
wards these effects are augmented up to 5 = 2 A, but 
for wider interfaces the self-polarization energy decays 
again. This singular feature can be understood if we re- 
member that the self-polarization energies actually are 
the product (integrated over the whole space) between 
the self-polarization potential V s and the electron (hole) 

probability density given by \ip e \ 2 f|V>fc| 2 ) [Eq. (|l|)]. 

To help with the explanation we present the Fig. 5 
where \ip e \ and V s are plotted separately. As can be im- 
mediately grasped from this figure, the maximum value 
for the self-polarization energy correction for the combi- 
nation of total confinement and step-like dielectric inter- 
face is a consequence of the fact that the wave-function 
is different from zero precisely in the same region where 
the self-polarization potential is positive (inside the dot). 
The situation changes drastically by relaxing the hard- 
wall boundary condition, as then | -0 e | 2 spreads consider- 
ably outside the dot, where the self-polarization potential 
is negative. This explains the strong decrease of E e + £/, 
by going from Vb e — * oo to Vb e = 4 eV; actually, as E^ 
is always positive (the holes are perfectly confined in the 
present case) an almost zero value of E e + £/j implies 
a negative value for E e . As the dielectric interface in- 
creases in size, the negative part of the self-polarization 
potential moves apart from the electronic interface, £ e 
becomes less negative, and E e + E/j increases. Finally, for 
<5 > 3 A, only the positive component of V s contributes 
to E e , and because its value decreases as 8 increases, E e 
shows a non-monotonic behavior. A final remark about 
Fig. 4 is that both Ecoui and E e + Eh values never exceed 
those achieved with the total confinement approximation, 
which seems to be a kind of upper limit. 
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FIG. 5. Upper pannel: probability density associated with 
the normalized ground state electronic wave-function, and 
effective mass mismatch mi n /rn out = 0.067. Full line, 
Voe — * oo; dashed-dotted line, Vbe = 4 eV. Lower pannel: 
self-polarization potential for 5 — 0, 0.5, 2, 5 A Dot radius, 
R= 10 A. 

Concerning the excitonic energy, the results seem to be 
less predictable. In Fig. 6 we show the exciton energies 
E ex (d) as a function of dot sizes, for quantum dots with 
S = 0.05, 0.5, 2, and 5 A; the IEMA results for S = 
are also plotted (thick full line). Contrary to Ecoui and 
E e + Eft., much higher excitonic energies than those re- 
sulting from IEMA can be reached for QD's with thin 
interfaces. The explanation for this is clear from the 
results presented in Fig. 4: while the self-polarization 
energy decreases strongly by going from Vb e oo to 
Vbe = 4 eV, the Coulomb energy remains essentially un- 
altered under such change in confinement. Accordingly, 
the exciton energy which is the difference between these 
two magnitudes, increases and becomes similar in magni- 
tude to the Coulomb contribution (roughly, E ex ~ Ecoui 
j 2). As the dielectric interface increases in size, the self- 
polarization energy first increases up to a critical <5, and 
then decreases: this non-monotonic behavior of E e + E/j 
is followed by a non-monotonic behavior for the exciton 
energy displayed in Fig. 6. It is interesting to note that 
for dielectric interface sizes of about a lattice parameter 
(2 A < 5 < 5 A), the excitonic energy correction is quite 
small for all d > 20 A. Under these conditions, a mea- 
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sure of the optical band-gap should be compared directly 
with the single-particle band gap, without excitonic cor- 
rections. 

It is worth noting the richness of behaviors than can 
be obtained for the dot size dependence of the exciton 
energy by the combination of partial confinement with 
a finite size for the dielectric interface. This should be 
kept in mind when correcting single-particle band gaps 
with the excitonic contribution, in order of be able of 
compare with optical experiments, as this is done usually 
by using the infinite barrier and step-like model for the 
dot-confined exciton (full thick line of Fig. 6). 



Figure 6 P.G. Bolcatlo and C.R. Proetlo 




FIG. 6. Exciton energies as a function of dot size for S = 
0.05, 0.5, 2, and 5 A. The full thick line corresponds to the 
infinite barrier limit for electrons (and 5 = 0), and the re- 
maining lines are for a barrier of 4 eV. As explained in the 
text, the curve with open points correspond to the exciton 
energy obtained by replacing Ei„ = 12.6 by Si n (d) as defined 
byEq. (20). 

Before conclude, let us take again the issue of the elec- 
trostatic screening in semiconductor quantum dots, dis- 
cussing it on the light of the preceding results. In the 
first place, we can include easily in our calculations a 
size-dependent dielectric constant for the dot. According 
to the generalized Penn model JlJ 



1 + (2a I d) 1 



(20) 



with ei n being the infinite dot size (bulk) dielectric con- 
stant, a = 2-kEf / E g k,F, and I = 2. E g is taken as the op- 
tical gap of the dot-acting semiconductor, and kp and Ep 
are determined from the density-of valence electrons no- 
For GaAs, taking E g = 5.2 eV,E3 and n a = 0.1774/A 3 , 
we obtain a — 8.01 A. Inserting this size-dependent di- 
electric constant in our calculations, we obtain the curve 
with open points of Fig. 6 (corresponding to S = 0.5 A). 
It can be seen that this size effect, although it exist, is 
quite small for dots with d > 10 A. 

The related issue of the incomplete screening at short 
distances in small semiconductor clusters can be dis- 



cussed along a similar line of reasoning. 
Eq. (20) as 



Generalizing 



1 



1 



t (r,d) 



1 



1 



r I a 



(21) 



and a is a screen- 
ed™ (oo) correspond- 



where e in (d) is given by Eq. (20), 
ing parameter. For d — > oo, and 
ing then to the bulk dielectric constant, the spatially de- 
pendent dielectpic function of Eq. (21) is that proposed 
by Hermanson.EJ The screening parameter a can be ob- 
tained from the criteria that the Fourier transform of Ea. 
(21) fits the dielectric function xif Walter and Cohen.Eil 
Following Oliveira and Falicovjij we took a ~ 0.58 A 
as the characteristic value for GaAs. Without any cal- 
culation, we can argue against the importance of this 
effect due to the following facts: i) the spatially depen- 
dent screening, as defined by Eq. (21) for a given dot size 
is extremely short-ranged, with a typical length scale of 
about an atomic Bohr radius, and ii) according to Eq. 
(20), Ei n (d — > 0) — > 1, and then the correction given by 
the second term in Eq. (21) is smaller for smaller dots. 
In other words, in this limit £i n {r, d — ► 0) — * £i n (d — * 0), 
and the incomplete short-range screening reduces to the 
size-dependent dielectric function. Let us estimate quan- 
titatively the importance of the second term in Eq. (21) 
on the ground state binding energy of unconfined (bulk) 
excitons: using first-order perturbation theory the cor- 
rection (in units of the bulk exciton binding energy Rq) 
is given by 



AE 



R 



8[£i 



Z (2 + a */ a y 



(22) 



with i?Q = 



e 2 /2s in aQ, aQ = h 2 s in /e 2 ^, and fx = 
m e mh/(m e + m/j). Inserting in Eq.(22) the GaAs pa- 
rameters of Table I, we obtain AE/Rq ~ 10~ 3 . Hav- 
ing discarded the importance of incomplete short-range 
screening for the bulk, let us estimate its relevance for 
quantum dots. To this end, we have calculated this cor- 
rection for the case of the ground state binding energy of 
an on-center donor impurity. This could be considered as 
the limiting case of m^/ma 3> 1, as in this case the exci- 
ton and the donor impurity problems are equivalent. In 
order to simplify matters, we take here the infinite barrier 
approximation. The results are given by the upper two 
curves in Fig. 4b. Of these two, the lower one (circles) 
corresponds to use in the calculation of EqouI (d) for the 
donor impurity only the first term in the right-hand side 
of Eq.(21), while the upper curve (squares) is the result 
by using Eq.(21) as it is. Once again, short-range effects 
are found quite small for dot sizes d > 20 A. 



IV. CONCLUSIONS 

In this work we have studied the dependence of 
Coulomb, self-polarization, and excitonic energies of 
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spherical semiconductor quantum dots with partial con- 
finement on dot size and in the presence of a smooth 
dielectric interface at the dot boundary. We have devel- 
oped a numerical method which allows the study of the 
electrostatic properties of finite-size dielectric interfaces 
of arbitrary shape; essentially, it consist in replacing the 
continuous profile of the dielectric interface by a series 
of discrete steps. In the limit of very large number of 
steps (several hundreds), our method provides the exact 
solution to the electrostatic problem of a graded three- 
dimensional dielectric interface. The motivation for this 
analysis is twofold: from the mathematical point of view, 
the combination finite barrier plus step-like model for the 
dielectric interface is ill-behaved, as the self-energy of the 
partially confined particle (electron, hole, or both) di- 
verges in this situation. From the physical point of view, 
the interface between two semiconductors (or between 
semiconductor and vacuum) is never perfectly sharp as 
the step-like model of the dielectric interface assumes. 
The finite size of the interface becomes then a crucial 
feature for a sensible calculation of the self-energy, al- 
though is not so relevant for the Coulomb energy. 

We have studied the excitonic energy for linear and 
cosine-like spatial profiles for the dielectric interface, 
and for different values of the interface width. Inter- 
estingly, the excitonic energy shows a pronounced and 
non-monotonic dependence on the interface width, be- 
coming quite small for intermediate dielectric interfaces, 
but becoming quite sizeable for narrow dielectric inter- 
faces. Our results highlights the importance of consider 
the finite size of the dielectric interface, as quite a va- 
riety of excitonic corrections can be obtained depending 
on the combination barrier height - interface width. 
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